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In this article are developed the results of my work (reference 1) 
"The Stability of Plates and Shells beyond the Elastic Limit." A 
significant improvement is found in the derivation of the relations 
between the stress factors and the strains resulting from the 
instability of plates and shells. In a strict analysis, the problem 
reduces to the solution of two simultaneous nonlinear partial differ- 
ential equations of the fourth order in the deflection and stress 
function, and in the approximate analysis to a single linear equa— ' 
tion of the Bryan type. Solutions are given for the special cases 
of a rectangular plate buckling into a cylindrical form, and of an 
arbitrarily shaped plate under uniform compression. These solutions 
indicate that the accuracy obtained by the approximate method is 
satisfactory. 


1. EXPRESSIONS FOR THE FORCES AND MOMENTS IN TERMS 
OF THE STRAINS IN THE MIDDLE SURFACE 


• On'a moveable Darboux trihedron, relative to which we shall 
study the element of the shell, we choose the xy plane to be 
tangent to the middle surface, and the x and y directions 
along orthogonal curves (fig, 1) . 

The state of stress of the element is determined by the 
tensor of the stress S. Its components Z_, Z . Z„, are 

small compared to X , Y , and X , that Is, each layer' of the 

<y tj 

shell element parallel to the middle surface is in a state of 
plane stress. The intensity of stress in this layer will be 


t *"Uprugo— plasticheskaya Ustoichivost Plasteen." Prikladnaya 
Matematika i Mekhanika X, igk6, pp 623-638, 
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The state of strain of the element is determined "by the 
components of the tensor of the strains e : TT , e yy-’ e xy’ 

since the shears e xz , ey Z are Small > hut the relation to the 
strain e zz may he found from the condition of constant volume 
of the element , 

e xx + e yy + e zz, ~ 0 (!*2) 

The intensity of strain in this layer of the material is given 
hy the formula 



2 

\/3 




\/ xx 


2 1 s' 

4. e e x © £ 

• ~yy + xx yy + k B xy 


(1.3) 


In agreement with the laws for the elasticity and plasticity 
of materials the- stresses and strains are connected hy the relations 
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a i 
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( 1 . 


Here Cj_ = a^(e^) is determined for each material as a function 
of e^. The properties of this function are as follows. Within 
the elastic .limit, that is, for o’ where 0 * is a physical 

constant, Hooke’s Law 0 ^ = le^ always holds. Beyond the elastic 
limit 0 ^ = <|> ( e-j^ ) is a certain curve (fig. 2). If at a certain 

instant of time there occur infinitely small variations from the 
state of .strain, that is, the quantities e . . , receive 

increments Se^ . , ., then the increments of stress helow the 
elastic limit are given hy the formulae (l,4) hy setting c i = Ee.^, 
Beyond the elastic limit the increments of stress, for > 0 are 

given hy formulae ( 1 . 4 ) in accordance with the curve = < *> (e^ ) , 
hut for Se-^ < 0 in*accordance with the law of unloading 0 ^ = Le^ . 
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The problem of the stability of shells (plates) is stated as 
follows: Given, a shell under a given system of applied forces 

and with the states of stress and strain known. Required, the 
critical value of the external forces for which, at the same time, 
there is equilibrium with other possible states of strain infinitely 
close to the original state. 

Let the change in the first and second quadratic form .of the 
middle surface of the shell relative to the given equilibrium 
position be characterized by the parameters e^, 2^, and 

y l' X 2* T where e.^, are length ratios and 2e^ is the 

shear in the middle surface in the x,y plane, and X-j , X g , r 

are changes in curvature and twist. According to the Klrchoff 
hypothesis the Increments in length and shear at a distance z 
from the middle surface will be 


8e 


xx 


S/ '-l 


Se 


yy 


ZX, 


Se 


xy 


= 2c. 


2 zx 


(1.5) 


We seek the stress increments corresponding to the strains (1.5). 
For this it is necessary to take the variations of relations (1.4), 
The variation of the intensity of strain may be found by making use 
of (1,5), but afterward, we write for the variation of the work of the 
internal forces 0^ 8c-. , in terms of the stress components. 


01 5e i - Zytexx + Y y 8o yy + X y 8e xy 


(1.6) 


We introduce the nondimens ional quantities 
■^•x -^y h 

X x * a — , y * = JL x * = JL x. * = - x. , x 
01 y y (T L 2 1 2 


!v t *-I t (i ’ t) 


where h is the shell thickness. Then in agreement with (1,1) 


*X* 2 + V 2 - V'V + 3X/ 2 = 1 
From (1,6) and (1,5) we have 

8e 1 =e - = 6 - zX 


( 1 . 8 ) 


k 


* 
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where 


e , X x » 5l + v *€ S + sy « 3 ** - vv * VV + a V r * 


x » X *x + t *X * ?V T 
x 1 y @ - v 


>(1.9) 


For the variation of formulae (1.4) we note that 



in which hy the properties of the curve » «(«*), — £ — > °* 

We denote hy z Q = z 0 * ^ the coordinate of the layer for which 

the intensity of strain is unchanged (Se^ = 0) during instability. 

It is clear that 


z 


o 




(l.io) 


The variations of formulae (1.4) have the form 



All quantities entering into the right-hand side of these equations 
except the strains and the curvatures are known, since the original 
state of stress of the shell (whose stability is sought) is supposed 
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given. The quantities E,. E’ •=» — , E MI = - — - are shown In figure 2 

©i dej^ 

as tangents of angles. Young’s modulus "being constant "but E* and 
E" depending on the state of stress. : ' 

Before instability the shell may find itself wholly "beyond the 
elastic limit, or it may have elastic regions, elasto~plastic regions, 
and purely plastic regions. If the state of stress is momentless , 
then the region of elasto— plastic strains, that is, the region where 
part of the shell thickness is elastic, part plastic, is absent. In 
this paper we confine ourselves to the detailed stability investigation 
of compressed plates in which the state of stress is always momentless 
before instability. Hence we shall suppose that in the shells 
considered below the region of elasto— plastic strain is missing, 
before buckling (this assumption is not essential). 

After instability tho region, of the shell where the stress was 
originally elastic will be, generally speaking, elastically deformed, 
since the strain variations are assumed infinitesimal. The region of 
purely plastic strain will be, generally speaking, resolved after 
buckling into two — one remaining purely plastic, the other elasto— 
plastic. Figure 3 shows a section normal to a shell with the three 
designated regions (tho plastic region after instability is shaded). 

Let the surface z - z a represent the boundary between the 

regions, one of which is elastic after instability, the other plastic. 

For its determination we shall suppose that in the elasto— plastic 
zone, the plastic zone adjoins the shell surface z = + ~ and the 

elastic zone which originates as a result of unloading adjoins the 
surface z =-|, 

In the region of elastic strain and in the zone of unloading 
( z ^ z Q ) formulae (l.ll) take the form 

SS X = E(s 1 -X lZ ) 5 S y = S(s 2 -X gZ ) 6X y = | E^ -Tz) ' (1.12) 

In the region of plastic strain, tod in the zone of loading ( g > 2 ^,) 
of the elasto— plastic region, these formulae may be presented in the 
fonrr * 


''"From (1.12) and (1.13) it is seen that the variations * SS^ ... 

on the boundary z Q are continuous in the case where the original state 
. of stress corresponds to the beginning of flow, and equally so when, as 
a result of variation, the state of stress changes in proportion to the 
original state (reference 2). 
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6S X = (E* -E") S x ^X( Z ~ z 0 )' + E*(e i: - X^) 

&S = (E« - E") S *X( z - z Q ) + - * z) I (1.13) 


8X y = (1* - E”) X y *X(z _ z Q ) + | S f ( s 3 - T z) 


We proceed to the derivation of expressions for forces and 
moments arising 'in the shell during instability. Eor their deter- 
mination we have 

, h ,nh ■ . 

2 / 2 


P?- 

«i - / h 53C X 42 8T 2 ■ / a ’ w y 42 83 * /_ h ■ SI y 42 
0 2 ^ 2 w 2 


h 

/•2~ 


^ h 

■ r? 


SMy = / h 8X r z dz SM 2 =/ h SY y z dz &H - j 

[I ~ 7 ? V o U 


/ 


h 

o 


SX z ■ dz 

a y 

p 


In the region of -purely plastic strains we obtain for the force in 
agreement with (1.13): 


i- (s Tl - | 8T 2 ) = 6l - %>8 x * e 


JL 

E’h 


7- (sw 2 - I Ml) - e 2 - X V s 

• h x ^ 


1 ss * I r_ - X»X *e (1.14) 

•J 


E T h 


D 

.J 


E’h 


and for the moments 

_4 

3D 

_4 

- ’ 3D 


;( 8 M i-| 5M 2 ) 

- (sMg - 1 r%) - 


:< x + x ’ s / x 


x 2 + x’s y *x 


— 6H 
3D 1 


T + x‘x *x (1.15) 
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1 

j 


vhere 


E*h 3 

9 


E* - E" 


(l.l6) 


In the r egion of. -purely elastic strains, formulae (1,1k) and 
(1.15) hold, only E« = E" = E, k* = 0. 

Thus In the two regions, the forces are linear functions only 
°~ e x-» € j-* an< ^ 2%, the middle surface shear, and the moments 

are linear functions only of the changes in curvature. 

In the region of elasto— plastic strain s, the stresses §X x . . 

have different expressions for z > z Q and for z £ z, Q . Hence, 

the Integrals in the expressions for the ^orces and moments must 
he split into two parts. For example, 

.h h 

r £ [)* a p 2 


5M ! ~ 2 6M 2 " /_ h BS x 


5S_ z dz 


h z dz + 


PS z dz 

3C 


in which for the region z Q ~>, z > — — we take SS^ according to 

j* 2 

(1.12) and for the region — ^ z > z Q , according to (1.13). As 
a result of these calculations we obtain for the forces 
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and for the moments 


12 

h 2 


- | SMg^ = E + E* + (E - E*)z o »3 X^* - |(E - E’)(l - z Q * 2 )e 


E* — E" 


+ C 1 - z o^) 2 (2 +" V') S x * X * 


i|(sM2-|sM^. 


1 + E f + (E - E*)z * 3 


X * 
2 


|(E - E')(l - V 2 )‘ 




+ — — (1 - z d *) 2 (2 + z Q *)S *X+ 


— 8H = — — 

h 2 3 


s 


+ E T + (E - E‘)z * 3 




(E - E*)(l - s 0 '' 2 )s, 


+ — — (1 - V >0 2 (2 + z Q *)X *X* 


(1 .18) 


The dependency between forces and strains is nonlinear, since ij 
z Q * enters into the formula and from (1.10) it depends on the strains.’ 
From this fact proceed all the difficulties of solution of problems 
in shell stability beyond the elastic limit. ,j 

Further, it is essential that the ordinate ' depending on 

both the changes in curvature X^, Xg, t and on the strains e^, 

6g, be expressed only in the changes in citrvature and the 

forces 5 Tq_, 5T 2 , 8S. Multiplying the first equation of (l.l?) 

by X x *, the second by Y *, the third by 3X * and adding, we 

*7 J 

get 


X(1 - z 0 *) 2 + kz Q * r- 4 


S X *5T 1 + S y *5T 2 + 3Xy*fiS 


0 


(1.19) 


EhX* 
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By introduction of the notation £ for the ratio of the thick- 
ness hp of the plastic layer to the thickness of the shell 


i - 


h P 


1 - z. 


( 1 . 20 ) 


and solving equation (I.I9) for £ , we get 

i 


E -yW (1 + cpj ^ 1 Z x)(l + 9) 

E - E" X 


( 1 . 21 ) 


where 


9 - 


X S**®]. + V OT 2 + 3 V 8S E - B" 

; , X = as 1 — 


da 3 

deT 


1 - X 


EhX# 


E 


E 


( 1 . 22 ) 


Formulae (1.17), (.1,18) are appreciably simplified (otherwise 
conserving the principal complications) . if we consider only the 
beginning of flow, that is, we suppose that the shell material 
before instability exceeds the elastic limit very slightly. In 
this case 


E’ = E X* X = 


E - E" 
E 


Therefore, in the notation of (1.20.) the corresponding formulae 
have the form for the forces ■ „ 


1_ 

Eh 


l“l - 1 bt 2 


) ‘ 


+ .~ S < 2 X 
2 x 

V + r V*®* 


Eh 




(1.23) 
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for the moments 

(m ± - | 8M 2 ^ = + xs x *c 2 (3 - 2 £ >x 

iL (8M a - | 8M^ - JC a + XS/? 5 (3 - 2t>4 

L. 5H .= - fl + XX *^ 2 (3 + 2£) x . • 

3D • 3 y 


(1.24) 


■where D is the usual stiffness for Poisson r s ratio equal to 1/2. 


2. THE STABILITY OF COMPRESSED PLATES 


Denoting the "bending of the plate during instability "by w(x,y) 
and the displacements of points in the middle surface projected in 
the x,y directions by u(x,y), v(x,y), respectively, we have 

expressions for the changes in curvature X^ X , t, and the 


strains e]_. 


e 2' 

X. = 

V 

5^w 

x = 

T ^W 

1 

Bx 2 

2 By 2 

Bx By : 

e, = 

Bu 



1 

Bx 

3 2 \By 





( 2 . 1 ) 


The forces applied in the middle surface before instability 
may be written in the following form: 


Ti - ^ lcr i y -x"' . 

and their projection on the 


Tg = hc t Y * S = hcrjXy* 

Z-axis after instability in the form 
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Therefore , the condition of equilibrium of all forces applied to an 
element and projected on the z-axis, gives 


^ 2 6M 1 - 

~W* 



3x $y 


d 2 8M - 

+ ha\.X 

Sy2 1 


0 


( 2 . 2 ) 


The condition of equilibrium of the middle surface forces after 
instability will be ■ 




5ST & 

~W 



( 2 . 3 ) 


Finally, the compatibility condition for the strains has the 

form 

c> 2 e 9 2 e 

i + ™_2 „ 2 1, o (2.4) 

Sy 2 cftc 2 dx dy 


The combination of differential equations (2.2), (2.3), and (2.4) 
is necessary and sufficient for the solution of the problem of 
stability, if the corresponding boundary conditions are set up.- 
Indeed, according to (1.14), or to (1.24) and (1.20), the strains 


: 1> 


may be expressed in terms of the forces 


8T 


1-’ 


5T, 


&S 


and the curvatures X (bending w), following which the moments 
SMp, 5Mg, SH are functions of these same four arguments. Thus 


the problem reduces itself to four differential equations with four 
unknown functions, of which (2.2) is of the Bryan type, a.nd ( 2 . 3 ), 
(2.4) are of the type of equations in plane problems. 


In the region of purely plastic strain of the plate, (that is, 
such that the whole thickness, plastic before instability, remains 
plastic after instability), the system of differential equations 
is resolved into two. For simplicity we consider only the case 
of the beginning of flow. Substitution of ifhe values of 5M^, 

8H from (1.15) into (2.2) gives a differential equation for w 
of the Bryan type: 



*v* + 2 


a 2 

ck <3y 


v + 




XX 


(2.5) 
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where, in agreement with (1.9) and (2.1), 


x. = x„* Sz. + gx * i 2 l. + y * ^ 


ox 2 


k y 


etc By By 


(2.6) 


The two boundary conditions on w agree with the usual 
boundary conditions for the Bryan equation.: 

Solving equations (1.14) for the strains, we get 

i=k( 6T i-i ST 2 

/ — /vyjLUJ.1 

\ xs * / 

S *8T-i + S *5T + 3S *5S 

x 1 / 9 J v 


= — jsT, - i 8T 0 ^ + 12 (s * 60 ?, + S, *8T + 3X *80 

(1 - X)Eh \ 1 5 7 ) 

XS * 


e = — (6T0 — J 3T-. ) + 

2 Eh \ 2 2 V (1 - X)Eh 


(2.7) 


„ 38S 

2 e = 2 — + 

3 Eh (1 


3XX * / \ 

I (S *OT, + s *8T 0 + 3Iy*8s) 

X)Eh V 3 " 7 *■ * J 


Equations ( 2 . 3 ) are satisfied if the stress function F is introduced: 


5T 1 B 2 F 


5T 


2 _ cFt 8s _ _ 5 s f 

Eh ^y2 Eh Bx 2 Eh ox ey . 


(2.8) 


following which, analogous to (2.6) we denote 


• ‘ B 2 F 

t = S v * 2 JL + S * ~ - 3 V — 

7 Bx 2 7 


8y r " 


Bx By 


(2.9) 


we obtain the compatibility condition- for strain in the .form ■ 

# . 

B 2 B 2 B 2 \ xt 

— s * + s * - 3 X * 

- 2 X Bx 2 7 Bx By ? / 1 — X 


V^F = 


8 y c 


(2.10) 
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In order to -write the boundary conditions f6r this equation, 
it is necessery : to compute the variations of the normal force 5T V , 
.-and of the tangential force- .' SSy : on;a certain. curvilinear contour 
in the middle surface of the opiate* • . : 


If the outward normal v and the tangent e to the contour 
constitute a ! coordinate system such that by rotation the positive 
direction of V coincides with that-of y and the positive . 
direction of- s- coincides- with that : of. x , and if the angle 
between the -normal and the x— axis is denoted "by a (fig. 4), . 
then our quantities have the known expressions 

' 8T., — ST: T - ; : 

5T V = + — — - cos 2a + 8S sin 2a 

s -.' 2 

v ■ > (2.11) 


- 5T p •• ' '■ 

8Sv = - sin 2a ~ SS cos 2a 

o 


The purely plastic region of the plate may be bounded by a 
contour., part of which coincides- with the ; boundary of the plate, 
the part adjoins the elasto— plastic region. For the formulation 
of the stability problem in the first part, the boundary conditions 
have the form ; 

5T V = SSy = 0 (2.12) 


and in the second part 8T V , 8S^ must he continuous. 

It- is easy to show that during -instability the entire .plate j 
may not-, remain in the purely plastic . state ; -that, is, an elasto— , 
plastic region may come into being. Indeed, going back, ve shall 
have the uniform boundary conditions (2.12) on all external edges 
of the plate. But the differential equations (2.3) and (2.4) for 
conditions (2.7) will be also... linear . and homogeneous and so will 
have the unique solution 




= 5S 




It follows from (2.7) that = e 3 

the basis of (1.9) and (1.10), z Q = 0. But 


= 0, from which on 
z = z Q is the 
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■boundary between the elastic and plastic zones through the thickness 
of the plate and the condition z Q = 0 specifies that the middle 

surface is this boundary. It follows that a given region of a plat© 
is not purely plastic, but e last o— plastic, which contradicts the 
assumption. 

During instability of a plate beyond the elastic limit it will 
either go completely over to the elasto— plastic state, or there will 
remain purely plastic regions in it, which are not diffused through- 
out the plate . 


In the region of elasto— plastic strains, equation (2.2) on the 
basis of expressions (1.24) may be presented in the form 



- X * + 2 
5x 2 * 


>2 

o x * + 2 

chc <V ^ - 5y2 


Y *1 X$ 2 (3 - 2£)X (2.13) 


in which, as in equations (2.5)* (2.10), the operator in parenthesis 
acts like a multiplier on the quantity to its right. 


The condition of compatibility of strain (2.4) on the basis of 
(1.23) has the form 


V^F = — 
2 



S * 




X 


(2.14) 


where the stress function F is determined by formulae (2.8). The 
value of £, the ratio of the thickness of the plastic layer to 
the plate thickness, enters into equations (2.13) and (2.14), 
therefore they show compatibility; this quantity £ is expressed 
by formula (1.21) in which the function cp is, if use is made of 
the notation (2.9) 


= 2 \ t 

b 1 - X x 


(2.15) 


Equations (2.I3), (2.14) agree with the corresponding equations 
(2.5) and (2.10) at the boundary of the purely plastic and the elasto- 
plastic regions. Indeed, at this boundary, besides continuity in 
the values of the forces 8Tp , SSy, the moments 8My , SHy ‘ 

(where SHy* is the rotational moment according to the boundary 
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conditions of Kirchoff), the bending w and the slope of the tangent 
plane, there must also hold the condition. 


or 


■ ; h p 

5=1 


(2.16) 


From (1.21) for this condition we have cp = — \ and t = — — Xh, 

2 . 

following which the remarked coincidence of the equations is easily 
shown. 

The. boundary conditions^ for equations (2.13), (2.1k) on the 
elasto— plastic part of the contour, coinciding with the plate contour, 
yield the usual requirement ST V = 5Sy = 0 and two conditions relating 
to the bending w. 

Condition (2.15) or 


t = - Xh (2.17) 


represents in itself the equation of the boundary between the purely 
plastic and the elastp— plastic regions. 

The possibility of purely plastic regions arising at the same 
with' "the elasto— plastic regions follows from the fact that, the 
value of 5 in agreement with (l.2l) and (2.15) may take on values 
not lying in the interval 1 >5 >0 . Certain examples are given 
below of exact solutions of the stability of plates., and, in 
particular, the problem of the compressed, plate freely supported 
along two sides; the edges of the plate near the free supports, 
after instability, remain in the purely plastic state. 


3. EXAMPLES OF EXACT SOLUTIONS OF PROBLEMS 
IN THE STABILITY OF PLATES 


The integration of the system of differential equations (2.13) 
and (2.1k) in the elasto— plastic region, and of (2.5) and (2.10) in 
the plastic region with an undetermined boundary between them given 
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by (2.16), is fraught with significant mathematical difficulties. 
As was shown in 1, the stability problem simplifies when the 
variations of the forces in the middle surface are zero every- 
where. In that case the relative thickness t of the plastic 
layer is a known function of the coordinates, since from (1.22) 
9=0 and consequently 


i = 



(3.1) 


If the state of stress of the plate before instability Is 
uniform, the value of £ will be. constant, since in (1.22) 

d^I 


Se'j' 


will be the same for the whole plate. 


We call those solutions of stability problems approximate 
for which the variations STt, 5T 2 , SS of the forces are 

identically zero. Thus, the equations (2.3) of equilibrium and 
the boundary conditions (2.12) are satisfied, but, except in 
special cases, the compatibility condition (2.4) Is not satisfied. 
The simplicity of such a solution arises from the fact that in 
equation ( 2 . 13 ) the value of t is known and given by formula 
(3.1), as a result of which this equation becomes linear with 
constant or variable coefficients. It closely resembles the 
equation for the elastic stability of an anisotropic plate. 


The exact solutions of the system (2.13), (2.14) are undoubtedly 
of interest in their own right, but for us they have significance 
because they can be made use of to estimate the degree of exactness 
of approximate solutions. 


We discuss a certain class of exact solutions of stability 
problems for uniformly compressed plates of arbitrary shape and the 
solution for a rectangular plate In the case when buckling into a 
cylindrical shape is possible. 


a. Stability of a Uniformly Compressed 
Plate of Arbitrary Shape (Fig. 4) 


In this case the state of stress of the plate before instability 
Ib uniform and given by the formulae 

X x = Y y = -°i’ X y = 0 (3,2) 


II 


ii ill mini 
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where is the compressive stress along the edge and is also the 

. uniform stress intensity at any point in the plate. The resulting 
stresses according to (1.7) and (l. 1 *-) will "be 

V-v -- 1 V = 0 s i* ■ v “ I (3 - 3) 

For the- values of X and t we have the expressions from (2.6) 
and (2.9) 

X = ~^w t » - i 9^F (3 .fc) 

2 


(3-5) 


class of exact solutions 


_ Xh ,, 2y 
~ 8~ C 


Equation (2.14) takes the form 

- Mi ^x) = 0 

Neglecting the harmonic function, we obtain a 


as a result of which the value of cp in (2.15) is expressed in terms 
of and from- (1.2) we find 



(3-7) 


The fundamental differential equation of stability (2.13) is now 
linear with constant coefficients and has the simple form 


1 - £ 2 (3 _ 2g) 


k 

V w + 



V 2 w = 0 


(3.8) 


Its solution has been much studied for different shapes of 
plates and for different boundary conditions, although in connection 
with the elastic stability of compressed plates. 
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The value of- f (3.7) Is. little different from. the approximation 
( 3 .1) , and characterizes the degree of; deviation : of the exact, solution 
from the approximate;, • • • 

In the general case we have from (3*5) 


t = 


\h 

8 


2 \ 
5 x + r x 


) 


(3-9) 


where F-, is an arbitrary harmonic function. For continuous circular 
plates, for example, . r is a constant. According to (2.15) anh (l.2l) 

we now have an expression for £ in terms of X 

(3-10) 



following which equation (2.13), having in the given case the form 

ha. 


5 ^ 


1 r(3 -25) 


X + — i X = 0 
D 


(3.11) 


has only one unknown function X, By use of relations (3-^) it W 
he integrated once 


i 1 

L 


31 

4 




, 2 ha . 

V w + — - w = 
D 


r 

2 


( 3 . 12 ) 


where is a new harmonic function, also a constant for continuous 

circular plates, insofar as w and V S w must he finite In the middle. 
Equation (3.12), in view of (3.10), may he solved for V w, after 
which the problem reduces to the integration of only one linear partial 
differential equation of the second order (for circular places) 


v' w = $ (w,r' 1 ,r 2 ) 

The stress function' F is now determined, in accordance with (3-9) 
and (3.4), from the Poisson differential equation 


* 
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U 2 x + r x ) (3.13) 


As we see, the problem of the stability of circular plates may 
be solved in comparatively simple fashion through to the end. The 
details of a similar calculation will be clarified below for the 
example of a rectangular plate compressed in one direction. 


b. Stability of a Rectangular Plate Under the 
Condition of Plane Strain (Fig. 9) 


Such a case occurs if, when a rectangular plate of length 2 is 
compressed in the x— direction, the width b in the y— direction cannot 
change as a result of walls along the boundaries y = 0 and y = b . 
The plane x = 0 shown in figure where C = 2c and L = 27, will 
evidently be a plane of symmetry of strains . 


We assume the buckling to result in a cylindrical shape. In such 
a case, according to the conditions of the problem, we have for the 
stresses before instability 


j 3 

7L, = -p Y six.' X v = 0 cr, = ~ P 
x ^2^ y a 2 


> (3.14) 


\/3 


X* = - ~ Y '* = - JL S * = ~ 3**0 

x \/3 Y \/3 x 2 -Y 


After buckling, w * w(x), e 2 = “ °* 

From equations (l .2h) we have 


SS * 0 OTp.o | ST 1 


Since, in accordance with the equations of equilibrium, ST-^ = const. > 
and ST-, = 0 from the condition at the edge x = ~, then we have 

the case = bT ? = 5S = 0. In consequence, the approximate 

solution, as was noted: at the beginning of 3, here becomes exact. 
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The thickness ratio £ for the plastic layer is a constant 
and is determined by formula (3.1). The stability equation (2.13) 
takes the form 


d fow + hp d 2 w _ q 

dx 4 J)\l - X5 2 (3 - 2g)] <*x 2 


If the relative Karman modulus, expressed by 

da. 


K = 



4(1 - X) 


(1 + V'l - X)‘ 


is introduced, then wo get from (3.1) 


X 


4(1 - \/k) j, \/k 
(2 - Vic) 2 ' " 2 


(3.15) 


(3.16) 


(3.17) 


following which we may simplify the expression for the parameter in 
equation (3*15) 

, £E = $£ (3.18) 

d[i - X? 2 (3 - 20 J Bk 

Since k = 1 up to the elastic limit, and k =.0 in a small 
area where there is flow of the material, and since the character- 
istic value of the parameter 7 must be the same in elastic and in 
plastic problems, then it follows from (3.18) that the critical 
stress corresponding to the small area of flow, is zero. 

It is interesting to note that the Karman problem may be 
considered as a limiting case of the stability of a rectangular 
plate compressed in one direction, of small width b, for which 
the parameter 7 will have the expression 


2 4hp 

y as 

3Bk 
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and consequently the critical stress is zero at the small area of 
flow.- As seen from the preceding and following examples of exact 
solutions, the total loss of load-carrying ability of a plate, 
predicted in the Karman problem, does not occur, generally speaking. 
This circumstance has already been noted (reference 1) . 


o. The Stability of a Rectangular Plate Compressed 
In One Direction (Fig. 5) 

We shall suppose that the rectangular plate, sufficiently long 
in the y-direction and compressed only in the x-directlon, buckles 
into a cylindrical shape. In this case 


X x = ^i Y y = X y ” 0 


(3.19) 


S * = -1 s ■* = - 

x y 2 


X * a -1 Y * = X * = 0 
x 7 7 


By the conditions of the problem, all sections of the plate 
y = const . remain plane after buckling and so we have 


e ^ - 9 e p " cons ‘k ■ 


( 3 . 20 ) 


on the basis of which from (1.24), SS = 0. Besides this, = 0 

from the boundary condition at the edges x = t — and consequently 
it follows from (2.3) that ST^ = 0 everywhere.^ 


Since there are no forces in' the y-direction, 
condition 



we must use the 


(3.21) 


From the second equation of the system (2.14-) we have 
dx-2 Eh 2 4 1 


(3.22) 
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since- X.= -X^,, It is not difficult to convince one ' s self , that 

(3.22) is the integral of equation (2*14) . The function cp , by 
which is found the value of £ from- (1.21)., here has the form 


1ST, 




cp = - 


* 2 £ 2 


(1 - l)!! 2 *! (1 - X)hX- 4(1 - X) 


(3.23) 


The bending moment in any section is 


SMq = -D 


l -f xg 2 < 3 - 2O 

4 


X 


(3.24) 


and so the boundary condition on the edges x = t — is X n = 0. 

2 


It is clear from (3-23) that, X^ cannot be zero in the elasto— 
plastic region since e p j- 0 (this follows from the constancy of 
sign of £ Xq, positive along the entire plate , necessitating 
e p / 0 to satisfy condition (3-22))., Thus the, elasto-plastic 
region does not go up to the edges of the plate and stops at the 

section x = — — , The region adjoining this to the edge will be 
2 2 

purely plastic. Indeed, since t X is positive, then is also 

positive. It follows from (2.7) that in the purely plastic and in 
the purely elastic regions the force ST p has the same sign as 

egj that Is, is a tensile force.. But if, to the plate, compressed.' 

beyond, the elastic limit in the' x~direction, there is applied a 
tensile force in the y-direction, then the plate remains in the 
plastic state. One may convince one's self of this by formally 
calculating the value of Seq according to (1,8), which at the 
edges is equal to eq, but the strain according to (2.7) 

is negative, and so the value of Seq will be positive, that is, 
plastic strains before buckling remain plastic after buckling. 


From (1.21) and (3.23) we now have- ■ 

hX n , _ 

— - = — , p(£) = - 4 + 8 £ - 3X£ 2 (3-25) 

^2 p(5) 


J 
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From this we find the lower limit to the value of £ (p > 0) 

V /. I. 3X. 




(3.26) 


The fundamental differential equation of stability (2.13) takes the 
form 


d2 

A:: 2 


1 ~ \ ^ 2 (3 - 2 ?) 


ha., 

i + - 




0 


(3.27) 


By introduction of the notation. 


Q(£) = 4 - 91 ? 2 + 6l£ 3 | = 


we write equation (3.27) in the form 


d2 a + h! . o 


2x 


c 

a ~ T 


d| 


2 P P 


(3.28) 


(3-29) 


where p is the basic parameter determining the critical stress 

p 2 = (3.30) 


D 


The integral of equation (3-29) may he obtained by quadratures. 
Through introduction of the notation 


pfH - o 4 - 12g •(. 12^ 2 - 3lg 3 

(h - 8£ + 3X? 2 ) 2 


(3.31) 


we obtain as a result 





(3-32) 
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In the purely plastic region we hare for the force 8Tg and the 
moment SM-^ in agreement with the results of 2 and (3.19): 



_ h(l - X) 
h ~ 3\ 



(3.33) 


The fundamental differential equation takes the form 



d£ 2 1+ - 3k 



(3.34) 


The solution, satisfying the condition X-^ = 0 at the end 1 = 1, 
is written in the form 


X, 


= C^e 2 sin 


n(i - P 
\A - 3X 


(3.35) 


in which as a result of symmetry we consider only deflections in the 
right half of the plate (x > 0) . 

For determination of the five undetermined constants namely, 
the three integration constants C^, C 2 , C^, the boundary 

coordinate a and the critical number |i. we may, besides 
equation (3,21), write four more, conditions: Conditions of symmetry 


|| = 0 (3.36) 

conditions at the boundary region 

I = a i = 1 (3.37) 


two continuity conditions, of moment and shear force, which in 
accordance with (3.25) and. (3*35) take the form 

0, Sin -p_ cos ikmgl . Shii (§) (3.38) 

3 N< /k _ 3X h(4 - 3X) n/ 4 - 3X V ft - 3X hP d (l) \ a£ /a 
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The constant e is not necessary and does not enter the 

c i : ; * . • V 

grn A.- 

conditions insofar as they are , independent of — and — . 

"" ' .-*2 • e 2 

By making use of the prescribed conditions and Introducing 
a new unknown the relative thickness of the plastic : layer 

at x = 0, we get for the values of |i and 1, — a (the relative 
length of the purely plastic part) the following formulae 


4M . n _ 1 4 L 
u ~ ■ # ■ 1 - a = ■ ■ ' « 

v/2(l - X.) 4 M 


(3.39) 


where L and M are the integrals 


I R ^o> x — 2? + X£ 2 /»B(5 0 5 .rp 

L = / dR, M=/ J 1 U . ffL - to h0 ) 

; / R ( 1) n/b(£ 0 ) ~ B (S) J B(l) v/B(? 0 ) - B(5) , 


in which the value of is determined by the relation 


cot 2 l — l ; \ ts 2(4 - 31) 
\/ 2 (l - \) 


E(r)-E(l)l (3; 41) 


As was already established, the value of 1 — a. is positive, 
therefore the integral L must be positive, and for this it is 
required that 1 - 2£ + \? 2 > 0, that is, - 




,/T 


•: (3.42) 


By considering the estimate (3.26), which is also reasonable for 
£ 0 , we' see that this quantity is contained within narrow limits and 
close to the approximate value' ('3.1). It follows from this that the 
critical stress will differ only slightly from.. the approximate vglue. 


2 6 


MCA TM No.' 1188 


d. Approximate Solution of the Problem for a Plate 
in a Uniform State of Stress Before Buckling 

In this case the stress components X x , Y y , and X y and the 
stress intensity are constant everywhere; the quantity \ 

will also be constant, and hence £ by (3.1)* 


The x btiiI y axes in a given case may be so chosen that 
the X stress is zero (principal axes of. stress). The fundamental 

y 

stability equation (2.13) takes the form 




i - |(1 - *)V 2 £ 


+ 2 


|d - «W 


o\r 

3x~oy“ 


1 - |(1 - k)Y. 


jc2 


y 


o w 


ji7 


ha. 

~ 


X * ^2 + Y * — 
X Sx 2 y c>y 2 _ 


(3.43) 


in which the generalized Karman modulus is introduced in accordance 
with formulae (3 .16) and (3.17), since the relation 


kC 2 (3 - 2t) - 1 - k 


holds . 

The coefficients in equation (3.43) are all positive, since 
the largest value of each of the quantities X^*, Y y * is 

2 

— — and 1 > k > 0. 

\/ 3 

Hence, the problem may be solved as a linear differential 
equation of the Bryan type with constant coefficients, and in 
difficulty is little different from the corresponding elastic 
case . 


Translated by E. Z. Stowell 
National Advisory Committee 
for Aeronautics 



NACA TM Ho. 1188 


27 


REFERENCES 


1. Ilyushin^ A. A. : TJstoichivost Plastinok i Obolochek za Predelom 

Uprugosti. Prikladnaya Maternatika i Melchanika, N. S. 8, No. 5 3 
19**, PP- 337-360. (Also available aa NAGA TM No. Ill6.) 

2. Ilyftshin, A. A.: K Teoria Malikh Uprugo— plastioheskikh Deformat sii. 

Prikladnaya Maternatika i Mekhanika, X, 19*+6, p. 3^7 • 








30 



Figure 5. 




